####################################################################################
###
###   Immigration, Public Housing and Support for the French National Front
###   Gloria Gennaro
###
###   Paper Figure A31
###
####################################################################################

library(rddensity)

################################################################################
# Set Up
################################################################################

# Load Working Directories
source(paste0(wd_main, '/2_code/00_working_directories.R'))

# Load data and clean
source(paste0(wd_code, '/05_data_load_clean_absun.R'))


################################################################################
# Plot Functions
################################################################################

denstest <- rddensity(df_did$running, p=2, c = 0)
rdplotdensity(denstest, df_did$running)
ggsave(paste0(wd_res, '/figures/figA31a.pdf'))

temp = df_did[abs(df_did$running)<2500,]
denstest <- rddensity(temp$running, p=2, c = 0)
rdplotdensity(denstest, temp$running)
ggsave(paste0(wd_res, '/figures/figA31b.pdf'))

temp = df_did[abs(df_did$running)<1500,]
denstest <- rddensity(temp$running, p=2, c = 0)
rdplotdensity(denstest, temp$running)
ggsave(paste0(wd_res, '/figures/figA31c.pdf'))

temp = df_did[abs(df_did$running)<1000,]
denstest <- rddensity(temp$running, p=2, c = 0)
rdplotdensity(denstest, temp$running)
ggsave(paste0(wd_res, '/figures/figA31d.pdf'))



